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I. INTRODUCTION. 



The structure of lamellar phases of symmetric AB diblock copolymers in a thin film is inves- 
tigated. We quantitatively compare the composition profiles and profiles of individual segments 
in self-consistent field calculations with Monte Carlo simulations in the bond fiuctuation model for 
chain length A'' = 32 and X-A = 30. Three film thicknesses are investigated, corresponding to parallel 
oriented lamellae with 2 and 4 interfaces and a perpendicular oriented morphology. Taking account 
of capillary waves, we find good quantitative agreement between the Monte Carlo simulations and 
,J^ ■ the self-consistent field calculations. However, the fluctuations of the local interfacial position are 

^ ' strongly suppressed by confinement and mutual interactions between lamellae. 
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'^ , Confinement of spatially structured fluids into thin films gives rise to a rich and interesting interplay between the 

intrinsic length scale of their structure in the bulk and the geometry of the film. One the one hand, confinement 

generally modifies the phase diagram of fluids m . For spatially structured phases the conflnement leads to transitions 

between phases with identical symmetry but different orientation with respect to the conflning walls. On the other 

i-rt ' hand the surfaces alter the local structure of the fluid in their vicinity. This gives rise to the enrichment of one 

r^ ■ component at the surface and the formation of wetting layers, but it also entails changes in the local fluid structure 

O ' (e.g., packing and alignment effects). Controlling the properties of confined fluids is important for various practical 

, O, applications. The orientation of the morphology is, for instance, important for connecting the self-assembled material 

to other devices (e.g., to form a molecular sieve 0). 

T— I : Studying effects of confinement on spatially structured polymeric fluids is particularly rewarding, because the large 

^ ' length scale of the phenomena, set by the chain extension, leads to a rather universal behavior, i.e., many features 

CO ■ are independent from the details of the chemical architecture. The equilibrium properties of polymeric systems are 

^~~^ ■ usually well describable by self-consistent field calculations. Moreover, the large length scales facilitate the application 

, -«^ ■ of experimental techniques. The detailed composition profiles in the vicinity of surfaces and even profiles of individual 

r-^ ■ segments are experimentally accessible. However, the large lengths scale of the phenomena goes along with protracted 

^\^ ' long time scales to reach equilibrium. There are quantitative comparisons between experiments and theory: The 

^^ , profiles of individual segments in films of diblock copolymers and their blends have been measured by small angle 

"^ ' neutron scattering, and the results have been compared to self-consistent field calculations y,g. The studies found 

Cd that the self-consistent field calculations describe the qualitative features, but that the profiles were broadened by 

H long wavelength fluctuations of the local position of the interfaces. 

^ ' In the previous paper Q, hereafter denoted as paper I, we have used a self-consistent field technique similar to 

Matscn Q for calculating the phase diagram of a symmetric AB diblock copolymer melt confined into a thin film as a 
function of the film thickness and temperature. The calculations have been compared to Monte Carlo simulations in 
^ , the framework of the bond fluctuation model at intermediate segregation {xN — 30) . Depending on the film thickness 
the systems assembled into parallel oriented lamellar phases L2 and L4, in which 2 or 4 AB interfaces run parallel 
to the surfaces of the film, or a perpendicular lamellar phase L± in which the interfaces run perpendicular to the 
confining surfaces. The dependence on the film thickness exhibited by the Monte Carlo simulation is in agreement 
5— j ' with the self-consistent field calculations. 

The aim of the present paper is to investigate the local structure of these films in more detail. The detailed profiles 
of the composition and individual segments provide a sensitive testing bed for comparing Monte Carlo simulations 
and self-consistent field theory. While the comparisons between experiments and theory have provided many valuable 
insights, Monte Carlo simulations might contribute to our understanding by investigating model systems in which many 
different quantities are simultaneously accessible. In the simulations we can analyze profiles of different quantities 
on various lateral length scales, and we can compare the AB diblock copolymer melt with a blend of A and B 
homopolymers of identical architecture. 

There are several potential sources of deviations between the Monte Carlo simulations/experiments and the self- 
consistent field calculations: (i) At low incompatibility, i.e., around the onset of spatial ordering, fluctuations become 



important. These shift the transition temperature from its mean field value x^ = 10-5 to higher values of the 
incompatibility, and the transition becomes first order ||7[] . Shifts of the transition temperature of similar magnitude 
are predicted in the P-RISM framework by Schweizer and co-workers Q . In our Monte Carlo simulations we observe 
the order-disorder transition in the regime 13.5 < x^ 5: 17 (cf. Fig. 7 in paper I). At higher segregation, however, 
fluctuation effects become less important, (ii) The self-consistent field theory assumes the interface between the A 
and B domains to be ideally flat in the lamellar phase. However, there are long wavelength fluctuations of the local 
interfacial position (i.e., capillary waves). Their free energy cost vanishes as their wavelength diverges, and, hence, they 
are important at all temperatures. These fluctuations lead to a broadening of the laterally averaged profiles measured 
in experiments or Monte Carlo simulations, (iii) Moreover, computer simulations P] and experiments |10-12J reveal 



that the copolymers assume a dumbbell-like shape already above the order-disorder transition temperature, as to avoid 
energetically unfavorable contacts between their different moieties. This leads to an increase of the chain extension. 
The interplay between the inter- and intramolecular energy has been studied in the P-RISM framework by Schweizer 
and co-workers [|j and field theoretical approaches [^ (iv) At high incompatibility (x ~ Oil)), the smallest length 
scale of the spatial structure (e.g., the width of the interface between the A and B domains) decreases rapidly. If this 
length scale becomes comparable to the length scale of the underlying microscopic structure (e.g., persistence length 
of the polymer or the range of packing effects of the monomers) a description in the framework of the Gaussian chain 
model breaks down. In fact, even for simple systems like interfaces between polymers of different persistence length 
the Gaussian chain model gives qualitatively erroneous predictions at high incompatibility |14| . Density functional or 
P-RISM calculations of more realistic models might account for the changes in the local fluid structure. Comparing 
self-consistent field calculations in the framework of the Gaussian chain model and density functional theory for a 
model with local fluid structure, Nath and co-workers ||l5|,|l6[ observed rather pronounced deviations between both 
approaches, (v) The confining surfaces create a sharp density gradient. In this region the local structure of the 
polymeric fluid is important at all incompatibilities. The monomer density exhibits pronounced oscillations in the 
vicinity of the walls (cf. Fig. 8 in paper I). Moreover, polymers align parallel to the wall. The latter effect is partially 
captured in the self-consistent field calculations and favors a perpendicular arrangement of the lamellae ||r^,|l3l . The 
aim of our paper to present a quantitative comparison between the Monte Carlo simulations and self-consistent field 
calculations and to quantify the effects discussed above in the framework of a well-studied model. 

The outline of our paper is as follows: In the next section we give a brief synopsis of the model and the computational 
technique. The reader is refered to paper I for further technical details. Then we discuss the strength of interfacial 
fluctuations and compare the laterally averaged profiles of the Monte Carlo simulations with those of the self consistent 
field calculations. Taking account of interfacial fluctuations we find good agreement, but the strength of the fiuctuation 
effect is smaller than expected. We close with a discussion of our findings. 

II. MODEL AND TECHNIQUES 

In the following we consider a melt of symmetric AB diblock copolymers confined into a thin fllm. The chain length 
is denoted by N . One half of the diblock consists of A monomers, while the other half consists B monomers. There is 
a short range repulsion x between A and B monomers, which drives the microphase separation. The lateral extension 
of the fllm is denoted as i, while Aq is the distance between the parallel, hard conflning surfaces. The monomer 
density in the middle of the film is denoted by p. 

We employ the bond fiuctuation model (BFM) for the Monte Carlo (MC) simulations. In this coarse grained lattice 
model an effective monomer blocks all eight corners of a unit cube from further occupancy. Wc work at a monomer 
number density of p = 1/16; a value which corresponds to a concentrated solution or melt. Monomers along the chain 
are connected via bond vectors from the set [2,0,0], [2,1,0], [2,1,1], [2,2,1], [3,0,0], [3,1,0] including all permutations 
and sign combinations. The large number of bond vectors allows 87 distinct bond angles and gives a rather good 
approximation of continuous space properties. Monomers interact via a square well potential which is extended over 
the nearest 54 lattice sites, which constitute the first neighbor shell in the monomer density pair correlation function. 
Monomers of the same species attract each other with strength — e while a contact in the range of the square well 
potential between unlike species increases the energy by an amount e. The surfaces are parallel and impenetrable. 
An A monomer in the two layers adjacent to the surfaces decreases the energy by e^, while a B monomer in this 
region increases the energy by the same amount. In the MC simulations we use chain length N ~ 32, and energy 
parameters e = 0.1769fcBr and e^, = O.lfcsT. This corresponds to intermediate segregation x^ — 30. The monomer 
wall interaction is rather weak; the A component does not wet the surface. The configurations of the polymers are 
updated via local hopping attempts of individual monomers, slithering snake-like movements and exchanges of the 
identity A ^ B oi the two blocks. This allows an efficient relaxation of the chain conformation. One Monte Carlo 
step consists of 3 slithering snake attempts per chain, 1 local hopping attempt per monomer, and 1 A ^ B flip per 



diblock. 

The self consistent field (SCF) calculations employ the Gaussian chain model. In order to map the bond fluctuation 
model onto the Gaussian chain model we choose the statistical segment length b such that R^ — 6^(iV— 1) w 289 equals 
the end-to-end distance of the chains in the Monte Carlo (MC) simulations in the athermal limit (e ~ 0). All distances 
are measured in units of the lattice spacing. In accord with previous studies we obtain b — 3.05 p9[ |. Monomers of 
different species repel each other via a contact interaction of strength x- This Flory Huggins parameter x is related 
to the potential well depth in the MC simulations via x — ^Zce/ksT, where Zc denotes the effective coordination 
number |pO| . The latter quantity has been extracted from the intermolecular paircorrelation function to be Zc = 2.65. 
Hence, we estimate that e = 0.1769kBT corresponds to x^ = 30 in the Gaussian chain model. This correspondence 
between the parameters of the bond fluctuation model and the Gaussian chain model has proven useful for studying 
the bulk and interfacial properties of binary homopolymer blends pl| ] and ternary mixtures of two homopolymcrs 
and a diblock copolymer [^ . The effect of the confining surfaces is twofold. On the one hand the monomer density 
decays continuously from its value p in the center of the film to zero in a narrow region of width A„ — 0.15i?e close to 
the surface. On the other hand, A monomers in this surface region are attracted while B monomers are repelled. The 
strength of this monomer wall interaction is chosen such that the surface energy in the MC simulations coincides with 
the value of the surface free energy in the SCF calculations in the limit that one species covers the surface completely 
and packing effects are neglected (cf. paper I). Using a SCF technique of Matsen y] we have numerically solved the 
model in mean field approximation. Further details of the MC simulations and the SCF calculations can be found in 
paper I. 

The phase diagram in the SCF calculations is presented in FigJ^. At high segregation, we observe an alternation of 
parallel oriented lamellae if the film thickness A matches a multiple of the bulk period Di,, and perpendicular aligned 
lamellae at intermediate values of the film thickness. This behavior corresponds to the parameters of the simulation 
XN = 30 and 1.71 < A/ Re < 3.2. At lower segregation or larger film thickness, the frustration of the parallel lamellar 
phase due to a mismatch between the intrinsic length scale of the lamellar order and the confinement decreases, and 
one encounters direct transitions between parallel lamellar phases. In this case the perpendicular lamellar phase L± 
and the adjacent parallel lamellar phase form a triple point. For very small film thickness A < 3-Dfc we find a second 
order transition between the perpendicular lamellar phase and the disordered phase. Note that there is no transition 
between the disordered state and the parallel ordered morphologies. Upon increasing the incompatibility the order 
propagates gradually from the surfaces into the bulk (cf. also Fig. 7 in paper I). 

III. RESULTS. 

Several independent systems are quenched from e = to e = 0.1769/c_bT. For film thickness Aq = 30 and Ap = 56 
they assemble into parallel oriented lamellar phases L2 and L4 with two or four AB interfaces, while we find a 
perpendicular oriented lamellar phase for film thickness Aq — 46. The Monte Carlo runs consists of at least 3.8 • 10^ 
Monte Carlo steps. We do not observe transitions between different morphologies within a simulation run and cannot 
rule out non-equilibrium effects completely. But the morphologies which are observed agree for all but one system of 
thickness Aq = 56 with the predictions of the self-consistent field (SCF) theory. Moreover, the structures obtained 
are free of defects and the laterally averaged composition profiles of independent quenches agree to a high accuracy. 
This is illustrated in Fig.y for Aq = 46, where 6 systems have assembled into a perpendicular lamellar phase with a 
repeat distance D — 1.936i?e. This value is about 6% larger than the prediction of the SCF theory. The quantitative 
agreement of the composition profiles of independent systems shows that each system has sampled the fluctuations of 
the internal interfaces appropriately. 

A. Interfacial fluctuations 

Fig.@ compares the two dimensional composition profiles in the perpendicular oriented lamellar phase at x^ — 30 
and A/JJe = 1.83. The walls (z-axis) run horizontally on the top and on the bottom of the figure. The x axis 
across the film is vertical. In each configuration of the MC simulations we have located the instantaneous position 
of the lamellae in the z direction, and shifted the z coordinate such that the center of an A lamella is located at 
z = 0. These shifted profiles have been averaged in the y direction and over all configurations. Hence, the profiles in 
the MC simulations are broadened by interfacial fiuctuations with wavevectors in the xy plane. The effect of these 
capillary waves is readily observed: the interfacial width in the MC simulations is significantly broader than in the 
SCF calculations. Similar deviations have also been observed in other studies [E3^,E4|. 



Apart from this broadening of the profiles, the MC simulations and the SCF calculations share many subtle details. 
In both cases the AB interface bends in the vicinity of the wall as to increase the surface fraction covered with the 
component of the lower surface free energy. The bending of the interface in the vicinity of the surfaces interferes with 
the distortion which originates from the opposite surface and leads to an oscillation of the interfacial position. This 
has been observed in SCF calculation ll^ and is also born out in the MC simulations. In both cases the interfacial 
width increases in the vicinity of the surfaces. In the SCF calculations this is due to a reduction of the density at 
the surface, while in the MC simulation the effect stems from the finite interaction range; monomers at the surface 
have less neighbors to interact with. This gives rise to a negative line tension when the AB interface approaches the 
surface ||]. 

In binary blends the interfacial structure can be characterized by the local interfacial position u{x,y), which depends 
on the lateral coordinates, and the "intrinsic" profiles of the order parameter across the interface. Neglecting the 
coupling between long wavelength fluctuations of the local interfacial position and the "intrinsic" profile, the latter can 
be calculated as a profile across an ideally flat interface. The SCF technique calculates these "intrinsic" profiles. The 
fluctuations of the local interfacial position in binary blends are well describable by the capillary wave Hamiltonian 

n[u{^)] = ^ld\[Vuf (1) 

where cFcS denotes the effective interfacial tension between the unmixed phases and u denotes the deviation of the local 
interfacial position from its lateral average. For unconfincd interfaces in blends the value of the effective interfacial 
tension agrees well with independent measurements of the free energy costs of an AB interface or SCF calculations 
|l4| . For interfaces in binary blends near a wall the effective interfacial tension extracted from the fluctuation spectrum 
depends on the distance between the interface and the wall ||2^, and the interaction between the interface and the 
wall imparts a long wavelength cut-off to the spectrum of interfacial fluctuations p^ . 

We investigate the fluctuations of the local interfacial position in more detail for the parallel lamellar phases L2 and 
£4. In the MC simulations we have measured the local interfacial position via a block analysis pq|. This is illustrated 
schematically in FigQa). We divide the simulation box into lateral columns of size B x B, and determine the local 
position of the interfaces in each column. The positions are Gaussian distributed around their lateral average across 
the system and s^ denotes the variance of the distribution. Using the capillary wave Hamiltonian one finds: 

2^Mi /L^ (2) 

where Xmax is the large length scale cut-off for interfacial fluctuations. If the interfaces were unconstrained this 
cut-off would be set by the lateral system extension Lmax = L [E8| . In the copolymer system neighboring interfaces 
interact and this might lead to a noticeable modification of the fluctuation spectrum. 

The results of this block analysis is presented in Fig.^(b) , which displays the variance of the interfacial position as 
a function of the lateral coarse graining size B. Circles refer to the L2 phase while squares correspond to the inner 
(in the middle of the film) and outer (close to the surfaces) interfaces of the L4 phase. The fluctuations of the outer 
interfaces in the ^4 phase and the interfaces in the L2 phase are similar. The spectrum of fluctuations is cut off at 
about Lniax ~ 40, indicating a rather strong interaction between the interfaces and the wall. The fluctuations for the 
inner interfaces of the L4 phase are well describable by Eq.(M). The spectrum is cut off by the lateral system size 
-^max — L = 96. From the slope of the straight line we estimate the effective interfacial tension to be: acS — O.lSTfcsT. 

This value can be compared to the interfacial tension of an AB interface in a binary blend [p9l aAB — jOfo-\/x/6(l — 
4 In 2/xN) « O-OGSiksT. Using a similar fluctuation analysis l^^l i the effective interfacial tension of a single unconfincd 
copolymer bilayer has been measured in the framework of the BFM at e = O.lSfc^T' or xN = 25. The analysis found 
that the effective interfacial tension cToff = O.OSfcsT of a copolymer bilayer with respect to undulations is much smaller 
than twice the interfacial tension in a binary blend 2aAB = 2 x O.OGGfcsT at that temperature. The latter finding 
is in accord with the expectation that the free energy cost of the interface in a copolymer melt is smaller than the 
interfacial tension in a blend. The interfacial tension in a blend is the limiting value at high segregation. 

Therefore, our finding for the confined copolymer melt suggests that interfacial fluctuations are strongly suppressed 
due to the confined geometry and mutual interactions between the interfaces. Estimating the additional free energy 
costs of interfacial fluctuations due to chain stretching, Semenov Q argued that equation (||) (with CTcff — <^ab) is 
appropriate when the spectrum is cut off at length scales larger than the distance between lamellae. Using (Jab, the 
width TTw of the AB interface as a lower cut-off and the periodicity D as an upper cut-off we obtain s = 0.16i?e. 
This value is larger than our MC results (cf. next section). Our findings suggest that the effective interfacial tension 
differs from the interfacial tension in a binary polymer blend. Other treatments of waves on the surface of brushes 
[pl|,p2[ result in wavevector dependent contributions to the free energy of fluctuations of the brush height. Laradji et 



al. |3^] have investigated fluctuations around the lamellar solution of the SCF theory in the bulk. Unlike the situation 
at an interface between partially miscible homopolymers, the eigenmodes of the fluctuations in the lamellar bulk 
phase are not capillary waves of individual interfaces but coherent modulations of the whole stack of interfaces. This 
points to a strong coupling between neighboring lamellae. Lacking a simple quantitative prediction for the fluctuation 
spectrum of an ensemble of strongly interacting interfaces in a confined geometry, however, we treat the strength of 
the interfacial fluctuations as an adjustable parameter in the following comparison. 

B. Profiles 

The comparison between the composition profiles in the L2 phase in the MC simulations and the SCF calculations 
is displayed in Fig.||(a). As discussed above, we find qualitative agreement, while there are quantitative deviations 
due to interfacial fiuctuations. To account for these effects, we assume that the intrinsic profiles averaged over a 
small lateral patch of the interface are describable by the SCF calculations while the local position of the interface is 
Gaussian distributed. To supplement the SCF results with this additional broadening we convolute the profiles p of 
any quantity with a Gaussian distribution [E3 34 : 



Pcupix) dx' pscf{x') exp{-[x ~ x'f /2s'^) (3) 

We employ the composition profile to adjust the width s of the distribution as to achieve best agreement between 
the MC simulations and the SCF calculations. We then use the same value of s to convolute profiles of all other 
quantities. Note that we assume the strength of the fiuctuations to be uniform across the film. Certainly, this is only 
a rough approximation. As we have observed in Fig.U(b) the fiuctuations of the inner lamellae are larger than the 
fiuctuations of those close to the film surfaces in the L4 phase. Hence, we expect s to be larger in the middle of the 
film and smaller at the surfaces. Therefore, by convoluting the SCF profiles we overestimate the fiuctuation effects 
at the film surfaces. 

Fig.^(b) shows the comparison between the MC simulations and the convoluted profiles of the SCF calculations 
with s/Re = 0.1. We achieve excellent agreement in the middle of the film and for the interfacial widths. The minor 
deviations at the surfaces are due to a reduced strength of interfacial fiuctuations and the different structure at the 
surfaces. The panels (c) and (d) present the comparison between the MC results and the unconvoluted and convoluted 
SCF profiles for individual segment densities. Again, the convolution improves the agreement between the MC results 
and the SCF calculations. The A ends are enriched at both surfaces, while the B ends are located mostly in the 
middle of the film. The maximum at the center is, however, not very pronounced. This indicates that the strong 
stretching limit is not yet reached for xTV = 30 and the two brushes formed by the copolymers largely interdigitate. 
The profiles of the A middle segment (monomer number 16) and the B middle segment (monomer number 17) exhibit 
a maximum at the interfaces. In accord with experimental observations, the distribution of the middle segments is 
narrower than the distribution of the chain ends [ p5[ . 

A similar analysis has been performed for the L4 phase and the L± phase. The results are presented in Figs.^ and 
0, respectively. The profiles in the L4 phase are qualitatively similar to those of the L2 phase. In the perpendicular 
phase L± the profiles run perpendicular to the direction of the AB interfaces of the lamellae and parallel to the walls. 
Since the A component is slightly enriched at both surfaces also in the B lamellae, the profile averaged across the 
film is less segregated in the B lamellae than in the A lamellae. This effect is observed in the MC simulations and 
the SCF profiles. 

For each phase we adjusted the strength of the interfacial fiuctuations as to match the SCF profiles onto the MC 
data. For the profiles of all other quantities we used the same value of s, and thereby improved the agreement between 
the MC simulations and the SCF calculations. This procedure results in s/Re = 0.1, s/R^ ~ 0.12, and s/R^ = 0.127 
for the 1^2 7^4 and L± phase, respectively. The increase of the strength of fluctuations is compatible with the intuition 
that the confined geometry reduces the interfacial fluctuations more in the L2 phase than in the L4 or L± phase. 
Using these values of s we can estimate the lateral block size Smin on which the MC data agree with the "intrinsic" 
profiles of the SCF calculations according to FigJ4(b). This yields the rough estimate Smin ~ 10. Semenov suggested 



this cut-off to be -Bmin — TTW for binary polymer blends and a recent MC study |34| in the framework of the BFM 
found Bmin = 1.27rw(l — 3.1/xN). Using w = b/^/6x, we obtain Bmin ~ 4.5. This is again an indication that the 
interfacial fluctuations are not well describable with the capillary wave Hamiltonian in the confined copolymer system. 
One possible explanation for a larger value of the cut-off Bmin is, e.g., a bending rigidity of the lamellae. 



IV. SUMMARY. 

We have presented SCF calculations and MC simulations for symmetric diblock copolymers confined into a thin 
film. Both surfaces attract the same component of the diblock via a short range potential. We have quenched several 
independent systems from the athermal state to xiV = 30, and we have compared the results of the MC simulations in 
the framework of the BFM with SCF calculations in the Gaussian chain model. We find qualitative agreement between 
the MC simulations and the SCF calculations. In particular, we observed the L2, Lj_ and L4 phases as predicted by 
the SCF calculations (cf. paper I). We have used these configurations to investigate the detailed structure in the thin 
film geometry. 

In the simulations we find evidence for a broadening of the profiles due to interfacial fiuctuations. However, the 
spectrum of interfacial fluctuations is not well describable by the capillary wave Hamiltonian. The effective interfacial 
tension is higher than the interfacial tension in a binary blend of A and B homopolymers, and the spectrum is cut — off 
at large wave lengths. A similar increase of the effective interfacial tension of an interface in a binary polymer blend 
in the vicinity of a wall has been observed [Em, however, the effect is more pronounced in the copolymer system. The 
confinement and the mutual interactions between neighboring lamellae strongly suppress interfacial fluctuations. 

To mimic the effect of interfacial fluctuations we convolute the SCF profiles with a Gaussian. The strength of 
the fluctuations is treated as a free parameter for each film thickness, however, we use the same value for profiles of 
different quantities. Taking account of interfacial fluctuations, we find almost quantitative agreement between the 
MC results and the SCF profiles for the composition, the density of chain ends and middle segments. In agreement 
with experiments |35| both MC simulations and SCF calculations show that the middle segements of the copolymer 
are stronger localized at the AB interface than the ends of the diblock in the middle of the domains. The good 
agreement between MC simulations and SCF calculations suggests that composition fiuctuations play only a minor 
role at intermediate segregation. Moreover, the local fiuid structure of the bond fiuctuation model does not have a 
large influence on the profiles. The latter finding is rather unexpected, because the interactions on the monomer scale 
X = 30/32 are large and the interfacial width between the A and B domains is not much larger than the microscopic 
length of the model. 
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FIG. 1. Phase diagram of a thin film with symmetric walls AiA'^ = A2A'' = 0.2 as a function of the incom- 
patibility xN and the film thickness A/ Re- L2, L4, and Le denote parallel lamellar phases with 2,4, and 6 AB 
interfaces, whereas L± denotes the perpendicular lamellar phase. The dashed lines mark multiples of the bulk 
lamellar period. The square denotes the approximate location of the triple point. 
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FIG. 2. yl density profiles of 6 independent quenches from x = to the L± phase at x^ = 30 for film thickness 
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FIG. 3. Contour plots of the composition in the perpendicular lamellar phase in a symmetric film A/R^ = 1.83 
and AiA'^ — A2N — 0.375 at x^ = 30 (a) SCF theory (lengths are measured in units of the end-to-end distance 
Re (b) MC simulations (lengths are measured in units of the lattice spacing/ Re — 17 lattice spacings). 
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FIG. 4. Fluctuations of the local interfacial position as a function of the lateral block length B: (a) Schematic 
illustration of the computational procedure for the L4 phase, (b) MC results for the parallel lamellar phases L2 
and I/4 in a symmetric film AiA'' = A2A'' — 0.375 at x^ ~ 30. 
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FIG. 5. Comparison between the results of the SCF theory and the MC simulations: L2 phase in a symmetric 
film /S./Re = 1.71 and AiA^ = A2N — 0.375 at x^ ~ 30. (a) composition profiles: symbols correspond to the MC 
results, lines represent the SCF calculations, (b) composition profiles: SCF profiles broadened with s — 0.0989. 
(c) segmental profiles, raw data: circles and triangles denote the middle monomers of the A block and the B 
block; diamonds and squares denote the end segments of the A and B block, respectively, (d) segmental profiles, 
SCF profiles broadened with s = 0.0989. Symbols as in (c). 
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FIG. 6. Comparison between the results of the SCF theory and the MC simulations: L4 phase in a symmetric 
film A/Re = 3.24 and AiN = A2A'^ = 0.375 at xN ~ 30. (a) composition profiles, symbols correspond to the MC 
results, lines represent the SCF calculations, (b) composition profiles, SCF profiles broadened with s — 0.12. (c) 
segmental profiles, raw data: circles and triangles denote the middle monomers of the A block and the B block; 
diamonds and squares denote the end segments of the A and B block, respectively, (d) segmental profiles, SCF 
profiles broadened with s = 0.12. Symbols as in (c). 
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FIG. 7. Comparison between the results of the SCF theory and the MC simulations: L± phase in a symmetric 
film A/Re = 1.83 and AiN = A2N = 0.375 at x^ ~ 30. (a) composition profiles, symbols correspond to the MC 
results, lines represent the SCF calculations, (b) composition profiles, SCF profiles broadened with s — 0.127. (c) 
segmental profiles, raw data: circles and triangles denote the middle monomers of the A block and the B block; 
diamonds and squares denote the end segments of the A and B block, respectively, (d) segmental profiles, SCF 
profiles broadened with s = 0.127. Symbols as in (c). 
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